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I— i. We present absolute parallaxes and relative proper motions for the central 

: stars of the planetary nebulae NGC 6853 (The Dumbbell) , NGC 7293 (The Helix) , 

^ . Abell 31, and DeHt 5. This paper details our reduction and analysis using DeHt 5 

Q_|. as an example. We obtain these planetary nebula nuclei (PNNi) parallaxes with 

Q ! astrometric data from Fine Guidance Sensors FGS IR and FGS 3, white-light 

^ ', interferometers on the Hubble Space Telescope (HST). Proper motions, spec- 

ci '. tral classifications and VJHKT2M and DD051 photometry of the stars com- 

! prising the astrometric reference frames provide spectrophotometric estimates of 

^ ! reference star absolute parallaxes. Introducing these into our model as obser- 

^ ! vations with error, we determine absolute parallaxes for each PNN. Weighted 

^ , averaging with previous independent parallax measurements yields an average 

! parallax precision, o-Tr/ir = 5%. Derived distances are: d-NGC 6853 = 405;1^25PC, 

g : d^GC7293 = 2 16l}^pc, dAbeU3i = 62ll?Jpc, and dDeHt5 = 345lJ?pc. These PNNi 

Q\ ', distances are all smaller than previously derived from spectroscopic analyses of 
the central stars. To obtain absolute magnitudes from these distances requires 
estimates of interstellar extinction. We average extinction measurements culled 



O 



^McDonald Observatory, University of Texas, Austin, TX 78712 

■^Centre for Astrophysics Research, STRI, University of Hertfordshire, CoUege Lane, Hatfield ALIO 9AB, 
UK 

^Department of Astronomy, New Mexico State University, Las Graces, New Mexico 88003 

^ United States Naval Observatory, Flagstaff Station, Flagstaff, AZ 86001 

^ Space Telescope Science Institute, 3700 San Martin Dr, Bahimore, MD 21218 

^ Department of Astronomy, University of Virginia, PO Box 3818, Charlottesville, VA 22903 

* Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 
16802 



-2- 



from the literature, from reddening based on PNNi intrinsic colors derived from 
model SEDs, and an assumption that each PNN experiences the same rate of 
extinction as a function of distance as do the reference stars nearest (in angular 
separation) to each central star. We also apply Lutz-Kelker bias corrections. The 
absolute magnitudes and effective temperatures permit estimates of PNNi radii, 
through both the Stefan-Boltzmann relation and Eddington fluxes. Comparing 
absolute magnitudes with post-AGB models provides mass estimates. Masses 
cluster around 0.57 M©, close to the peak of the white dwarf mass distribu- 
tion. Adding a few more PNNi with well-determined distances and masses, we 
compare all the PNNi with cooler white dwarfs of similar mass, and confirm, as 
expected, that PNNi have larger radii than white dwarfs that have reached their 
final cooling tracks. 

Subject headings: astrometry — interferometry — planetary nebulae — stars: 
distances — stars: white dwarf — stars: masses 



Introduction 



Planetary nebulae (PNe) are a visually spectacular and relatively short-lived step in 
the evolution fr o m asy mptotic giant branch (AGB) stars to a final white dwarf (WD) stage. 
Iben fc Renzinil (119831 ) first argued that the ejection of most of the gaseous envelope in AGB 
stars occurs during the thermal pulse phas e, in the form of a massive, low-velocity wind. 
As summarized by IStanghellini et al.l (120021 ). the remnant central star (PN nucleus, PNN; 
PN nuclei, PNNi) ionizes the gaseous ejecta, while a fast, low mass- loss rate PNN wind 
shapes the PN. PN morphology depends on a complicated combination of phenomena, some 
occurring within the nebular gas, which evolves on a dynamical timescale, and others caused 
by multiplicity and/or the evolution of the stellar progenitors and of the PNN. Morphology 
may also depend on the physical status of the interstellar environment of the PN progenitor. 

Intercomparison of PNe can aid our understanding of the complicated astrophysics of 
this stage of stellar evolution, particularly if distances are known. M any indirect metho ds 



of PN distance determination exist (e.g., ICiarduUo et al.l (119991 ) and iNapiwotzkil (120011 )) 



including estimates from interstellar Na D lines, NLTE stellar atmospheres analyses of the 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555 
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PNNi (e.g. Hultzsch et al. 2007 ), est imates from resolv ed companion stars fiCiardullo et al. 
19991 ). and from Galactic kinematics (iNapiwotzkil |2006| ) . The expansion method (e.g. Palen 
et al. 2002 ) becomes mo del- dependent when applied to PNe with asymmetric or irregular 
geometry, and requires either an assumption that apparent expansion is due to material 
motions, or that one models the motion of the ionization front. Agreement among these 
distance determination methods is seldom better than 20%, often worse. Direct parallax 
measurments of PNNi rarely have p recis ions smaller than th e measured parallax, a notable 
exception being [Harris et al.l (119971 ) and [Harris et al.l (120071 ). who, using narrow-field CCD 
astrometry, provide ~0.4 millisecond of arc (mas) precision parallaxes for PNNi nearer than 
-500 pc. 

To further reduce the distance errors for a few PNNi (chosen as nearby from Harris et 
al. 1997), we have determined new absolute parallaxes of the PNNi of DeHtS, Abell31, and 
NGC 7293 (The Helix), using FGSlr. Positions and aliases are given in Table 1. We have also 
determined a revised parallax for NGC 6853 (The Dumbell = M27), using previously collected 
FGS 3 data (Benedict et al. 2003), combined with several new FGSlr measurements. Our 
present errors average 2-3 times smaller than those in Harris et al (2007). However, some 
reduction in the final pa rallax errors i s obta in ed through a weighted average of o ur present 
parallaxes with those in [Harris et al.l (120071 ). iNapiwotzki fc Schoenbernerl (119951 ) classifies 
all the PNNi considered in this paper as WD of type DAO or DA. 

Our reduction and analysis of these data is basically the same for our previous work on 
NGC 6853 (Benedict et al. 2003). Our extensive investigation of the astrometric reference 
stars provides an independent estimation of the line of sight extinction as a function of dis- 
tance for these PNNi, a significant contributor to the uncertainty in the absolute magnitude. 
My. Using DeHt 5 as an example throughout, we present the results of extensive spectropho- 
tometry of the astrometric reference stars, information required to derive absolute parallaxes 
from relative measurements; briefly discuss data acquisition and analysis; extract limits on 
binarity and photometric variability; and derive an absolute parallax for each PNNi. Finally, 
from a weighted average of our new parallaxes and those of Harris et al. (2007) we calculate 
an absolute magnitude for each PNN and derive stellar radii. With these and estimates of 
PNN mass from po st-AGB evolution m odels we derive surface gravities, log to compare 
with those from the INapiwotzkil (119991 ) stellar atmosphere analyses. We discuss some astro- 
physical consequences of these new, more precise distances and summarize our findings in 
Section [71 



Bradlev et all Jl 9 9lh andlNelanI (j20o3) pr ovide an overview of t h e FGS mstrume r it and 



Benedict et al. I (Il999f ) 



Benedict et al. 



(l2002bl ). [Harrison et al.l (120041 ) . Benedict et all (l2007l ) 



describe the fringe tracking (POS) mode astrometric capabilities of an FGS, along with the 
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data acquisition and reduction strategies used in the present study. We time-tag all data 
with a modified Juhan Date, mJD = JD - 2400000.5. 



Observations and Data Reduction 



Using DeHt 5 as an example, Figure [T] shows the distribution on the sky of the ten 
reference stars and the PNN. This image is from the Digitized Sky Survey, via Aladin. For 
this target fifteen sets of astrometric data were acquired with FGS Ir on HST, spanning 3.8 
years, for a total of 248 measurements of the DeHt 5 PNN and reference stars. Each data set 
required appr oximately 33 mi n utes of s pacecraft time. The da t a were reduced and calib rated 
as detailed in [Benedict et all J2002bh . [Benedict et all J2002ah . [McArthur et all J200lh . and 
Benedict et air ( 2007 ). At each epoch we measured reference stars and the target multiple 
times, this to correc t for intra-orbit drift of the type seen in the cross filter calibration data 
shown in figure 1 of [Benedict et al.l (j2002bl ). 



Table [2] lists the fifteen epochs of observation and highlights another particular difficulty 
with these data. Ideally (cf. Benedict et al. 2007) we obtain observations at each of the two 
maximum parallax factorq^ at two distinct spacecraft roll values imposed by the requirement 
that HST loW to provide thermal control of the camera in the radial bay and to keep its solar 
panels fully illuminated throughout the year. This roll constraint generally imposes alternate 
orientations at each time of maximum positive or negative parallax factor over a typical two 
year campaign. A few observations at intermediate or low parallax factors usually allows a 
clean separation of parallax and proper motion signatures. In the case of DeHt 5, (as well as 
Abell31 and NGC7293), two-gyro guiding forced us into the less than satisfactory distribu- 
tion of parallax factors shown in Tabled Specifically for DeHt 5, there are no large positive 
parallax factors in Right Ascension. However, the higher the absolute declination of the 
target, the more likely it is that there will be windows of visibility near times of ± maximum 
parallax factor, either in RA or declination. Additionally, large declination typically results 
in higher ecliptic latitude. The ecliptic latitude of DeHt 5 renders its parallactic ellipse rather 
round, increasing the value of observations that were forced to be secured at times far from 
maximum parallax factor. We gain parallactic displacement in declination at the expense of 



^Parallax factors are projections along RA and Dec of the Earth's orbit about the barycenter of the Solar 
System, normalized to unity. 

^HST has a full compliment of six rate gyros, two per axis, that provide coarse pointing control. By 
the time these observations were in progress, three of the gyros had failed. HST can point with only two. 
To "bank" a gyro in anticipation of a future failure, NASA decided to go to two gyro pointing as standard 
operating procedure. 
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displacement in RA. These scheduling and solar-panel illumination constraints also resulted 
in the spotty access to our reference stars indicated in Table |2l 

Finally, for DeHtS and AbellSl we were able to take advantage of science instrument 
command and data handling (SIC&DH) computer problems that took the only other then 
operational science instrument (WFPC2) off-line in late 2008. This situation opened a 
floodgate of FGS proposals, temporarily rendering HST nearly an 'all astrometry, all the 
time' mission. Consequently, we obtained another epoch well-separated in time from the 
original eleven. This permitted a significantly better determination of relative proper motion 
for these two targets. 



3. PNNI Photometry and Companion Limits 

The FGS have two operating modes. We used both. The fringe tracking (POS mode) 
data collected over the course of this project can be used to establish the degree of photomet- 
ric variability of these PNNi. Additionally, fringe scanning (TRANS mode) data provide an 
opportunity to discover previously unknown companions, or to establish limits for separation 
and magnitude difference. 



3.1. FGS Photometry of the PNNi 

FGS Ir and FGS 3 are precision photometers, yielding relative photometry with internal 
error > 0.002 mag (Benedict et al. 1998a). During each of the observation sets we observed 
the PNNi 4-7 times over approximately 33 minutes. Given the faintness of these stars, we 
made no effort to explore for high-frequecy variations during a single observation]^ Standard 
deviations within any one observation sets were on order 0.1% for the two brighter PNNi 
(NGC6853 and NGC7293), 0.2% for Abell31, and 0.5% for DeHtS. We derived a PNNi 
average intensity for each observation set and used the summed intensity of the brightest 
astrometric reference stars as a flat field. The resulting intensities were then transformed 
to relative magnitudes such that the a verage relative mag nitude over the entire campaign 



matched the measured magnitude from [Harris et al.l (120071 ). 



We show a montage of the resulting PNNi light curves in Figure [2l Our coverage is 
too sparse to extract any believable periodic astrophysical component to these variations. 



*The FGS samples the fringe zero crossing at a 40Hz rate. See Benedict et al. ( 1998bh for an example of 



the use of an FGS as a high-speed photometer, capturing a flare event in the vicinity of Proxima Gen. 
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However, the internal errors for three PNNi suggest variations far larger than those expected 
from photometric errors only. This was expected. We fit these variations with sinusoids with 
periods of one year, because of previously encountered position dependent variations of the 
brightnesses of the stars used to generate the flat field (see Benedict et al. 1998a, section 2.3.1 
for a discussion of roll-induced photometric variation). All PNNi were observed at similar 
positions within the FGS, typically with radial distances from pickle center <10". Variations 
introduced by Zodiacal light will be insignificant for the two brightest PNNi and for DeHt 5 
at a high ecliptic latitude. AbellSl is closest to the ecliptic, but nearly all observations were 
obtained at essentially identical Sun-target angular separations, resulting in little variation 
due to changing background. From Figure [2] and the statistics of intra-orbit observations we 
conclude that an upper limit to photometric 'flickering' on timescales of minutes to years is 
around 5 mmag. 



3.2. Assessing PNNi Binarity 



Companions of stellar and substellar mas s have been i nvoked to produ ce the asymmet- 
ric structure seen in PN f lSond fc Liviolll99ol : ISokerl bood : IPe Marcol 120091 ) . Using FGS Ir 
TRANS mode observations (e.g., Franz et al. 1998 , Nelan et al. 2004 ), we have analyzed 
the fringe morphology of these PNNi and find that all are unresolved. This places limits on 
separation and magnitude difference. Am, for possible companions. Details on detectability 
can be found in Section 3.3.2 of the FGS Instrument Handbook (Nelan 2007). Summariz- 
ing the complex interplay among system magnitude. Am, and separation, we would have 
detected companions with separations of 10 mas and larger with Am < 1. Detectability 
at separations of 15 mas increases to Am < 2. For separations > 50 mas FGS Ir achieves 
a detectability of Am < 3.5. Figure [3] compares the X and Y axis fringes of DeHt 5 with 
those from Abell31. Either the two objects have identical companions, or they are both 
unresolved. Similar comparisons were made between NGC 6853 and NGC 7293 with similar 
results. Once we determine absolute magnitudes, parallaxes, and estimate masses for these 
PNNi we will (Section 16.5.31 below) establish the spectral types, separations in AU, and 
periods for companions that would remain hidden from the FGS. 



4. Spectrophotometric Absolute Parallaxes of the Astrometric Reference 

Stars 

Because the parallax determined for the PNNi will be measured with respect to reference 
frame stars which have their own parallaxes, we must either apply a statistically derived 
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correction from relative to absolute parallax (Van Altena, Lee & Hofileit 1995, hereafter 
YPC95 ) or estimate the absolute parallaxes of the reference frame stars listed in Table [2j 
In principle, the colors, spectral type, and luminosity class of a star can be used to estimate 
the absolute magnitude. My, and V-band absorption. Ay. The absolute parallax is then 
simply, 

■^Q-{V-Mv+5-Av)/5 



T^abs 



(1) 



The luminosity class is generally more difficult to estimate than the spectral type (tem- 
perature class). However, the derived absolute magnitudes are critically dependent on the 
luminosity class. As a consequence we obtained additional photometry in an attempt to con- 



f irm t he luminosity classes. Specifically, we employ the technique used by iMajewski et al. 
( I2OOOI ) to discri minate between giants an d dwarfs for stars later than ~ G5, an approach 
also discussed by lPaltoglou fc Belli (119941 ). 



4.1. Broadband Photometry 

Our band passes for reference star photometry include: BV (CCD photometry from a 
Im telescope at New Mexico State University) and JHK (from 2MAS^. We also had access 
to Washington/DDO filters M, T2, and DD051 (obtained at Fan Mountain Observatory 
with the Im, and at Las Campanas Observatory with the Swope Im). Table [3] lists the 
visible and infrared photometry for the DeHtS reference stars, ref-2 through ref-11. 



4.2. Spectroscopy, Luminosity Class-sensitive Photometry, and Reduced 

Proper Motion 

The spectra from which we estimated spectral type and luminosity class come from 
the New Mexico State University Apache Point Observatorj]^. The dispersion was 0.61 
A/pixel with wavelength coverage 4101 - 4905 A, yielding R~3700. Classifications used a 
combination of template matching and line ratios. The brightest targets had about 1500 
counts above sky per pixel, or S/N ~ 40, while the faintest targets had about 400 counts 
per pixel (S/N ~ 20). The spectral types for the higher S/N stars are within ±1 subclass. 



^The Two Micron All Sky Survey is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/California Institute of Technology 

^ The Apache Point Observatory 3.5 ni telescope is owned and operated by the Astrophysical Research 
Consortium. 
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Classifications for tlie lower S/N stars are ±2 subclasses. Table |3] lists the spectral types and 
luminosity classes for our reference stars. 

The Washington/ DDO photometry can provide a possible confirmation of the estimated 
luminosity class, depending on the spectral type and luminosity class of the star (later than 
G5 for dwarfs, later than GO for giants). Washington/DDO photometry was more helpful as 
a discriminator for the DeHtS field than for our previous work on NGC6853 (e.g., Benedict 
et al. 2003), suggesting a giant luminosity classification for ref-9 and ref-11. However, the 
problems related to the NGC6853 PN nebulosity discussed in that paper also affected the 
NGC 7293 reference frame DDO photometry. 

We employ the technique of reduced proper motions to provide a confirmation of the 
reference star estimated lu minosity class listed i n Table Ei We obtain preliminary proper 



motions (/i) from UCAC2 (jZacharias et al.ll2004l ) and J, K photometry from 2MASS for a 



one-degree-square field centered on DeHt 5. With final proper motions from our astrometric 
solution (Section 15.11) we plot Figure HJ which shows Hk = K + 51og(/i) versus (J — K) 
color index for 458 stars. If all stars had the same transverse velocities. Figure H] would be 
equivalent to an H-R diagram. DeHt 5 and reference stars are plotted as ID numbers from 
Table O Errors in Hk are now ~ 0.3 mag. Reference stars 9, and 11 are clearly separated 
from the others, supporting their classification as giants. 



4.3. Interstellar Extinction 



To determine interstellar extinction we first plot these stars on several color-color dia- 
grams. A comparison of the relationships between spectral type and intrinsic color against 
those we measured provides an estimate of reddening. Figure [5] contains a J-K vs V-K 
color-color diagram and reddening vector for Ay = 1.0. Also plotted are m appi ngs b e tween 
spectral type and luminosity class V and HI from iBessell fc BrettI (119881 ) and ICoxl (120001 ) 
(hereafter AQ2000). Figure and similar plots for the other measured colors, along with 
the estimated spectral types, provides an indication of the reddening for each reference star. 

Assuming an R = 3.1 galactic reddening law (Savage & Mathis 1979), we derive Ay 
value s by comparing t he me asured colors (Tabled) with intrinsic B-V, J-K, and V-K colors 
from iBessell fc BrettI ( 1l988l ) and AQ2000. Specifically we estimate Ay from three different 
ratios, each derived from the Savage & Mathis (1977) reddening law: Ay/E(J-K) = 5.8; 
Ay/E(V-K) = 1.1; and Ay/E(B-V) = 3.1. The resulting average Ay are collected in Table 

El 
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4.4. Adopted Reference Frame Absolute Parallaxes 

We derive absolute parallaxes with My values from AQ2000 and the (Ay) derived from 
the photometry. Our parallax values are listed in Table [31 We produce errors on the absolute 
parallaxes by combining contributions from uncertainties in My and Ay, which we have 
combined and set to 0.5 magnitude for each reference star. Individually, no reference star 
parallax is better determined than ^ = 23%. The average absolute parallax for the DeHt 5 
reference frame is (iTabs) = 1-0 mas. As a sanity check we compare this to the correction to 
absolute parallax discussed and presented in YPC95 (section 3.2, fig. 2). Entering YPC95, 
fig. 2, with the DeHt 5 galactic latitude, 1 = -12°, and average magnitude for the reference 
frame, {Vref)= 14.3, we obtain a galactic model-dependent correction to absolute of 1.0 mas, 
in agreement. 



5. Absolute Parallaxes of the PN Central Stars 

Sections 5.1-5.4 detail our astrometric modeling of the DeHt 5 data. Any differences in 
modeling for other PNNi are noted in Section 15.51 below. We compare our new distances 
with other more indirect estimates later in Section 16.5.41 



5.1. The DeHt 5 Astrometric Model 

With the positions measured by FGS Ir we determine the scale, rotation, and offset 
"plate constants" relative to an arbitrarily adopted constraint epoch (the so-called "master 
plate") for each observation set (the data acquired at each epoch). The mJD of each ob- 
servation set is listed in Table [2l The DeHt 5 reference frame contains 10 stars. We employ 
an eight parameter model for those observations. For the DeHt 5 field all the reference stars 
are redder than the science target. Hence, we also apply the corrections for lateral color 
discussed in Benedict et al. (1999). 



As for all our previous astrometric analyses, we employ GaussFit ( iJefferys et al.lll988l ) 
to minimize x^- The solved equations of condition for DeHt 5 are: 

x' = X + lc^{B - V) (2) 

y' = y + lCy{B-V) (3) 
^ = Ax' + By' + C + R,{x'^ + y'^) - fi^,At - P.vr, (4) 
r] = Dx' + Ey' + F + Ry{x'^ + y'^) - fiyAt - Pgny (5) 
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where x and y are the measured coordinates from HST; lex and Icy are the lateral color 
corrections from Benedict et al. 1999; and B — V are those colors for each star. A, B, D, 
and E are scale and rotation plate constants, C and F are offsets; and Ry are radial term 
coefficients; /i^ and /i^ are proper motions; At is the epoch difference from the mean epoch; 
Pa and Ps are parallax factors; and tt^ and vr^ are the parallaxes in x and y. We obtain the 
parallax factors (projections along RA and Dec of the Earth's orbit about the bary c enter of 
the Solar System normalized to unity) from a JPL Earth orbit predictor (IStandishI (Il990l )). 
upgraded to version DE405. 



5.2. Prior Knowledge and Modeling Constraints 

In a quasi-Bayesian approach the reference star spectrophotometric absolute parallaxes 
(Table [3]) and proper motion estimates for DeHt 5 (Harris et al. 2007) and for the reference 
stars from UCAC2 (Zacharias et al. 2004) were input as observations with associated errors, 
not as hardwired quantities known to infinite precision. Input proper motion values have 
typical errors of 4-6 mas y~^ for each coordinate. The lateral color calibration and the B-V 
color indices are also treated as observations with error. Orientation to the sky is obtained 
from ground-based astrometry from 2MASS with uncertainties in the field orientation ±0?05. 
This value, too, was made available as an observation with error. We essentially model a 
3^ volume of the space that contains our science target and reference stars, all at differing 
distances. 



5.3. Assessing Reference Frame Residuals 

The Optical Field Angle Distortion calibration (McArthur et al. 2002) reduces as-built 
HST telescope and FGS Ir distortions with amplitude ~ 1" to below 2 mas over much of the 
FGS Ir field of regard. From histograms of the reference star astrometric residuals (Figure [6]) 
we conclude that we have obtained satisfactory correction in the region available at all HST 
rolls. The resulting reference frame 'catalog' in ^ and r] standard coordinates (Table Hj) was 
determined with < cr^ >= 0.3 and < cr^ >= 0.3 mas. Relative proper motions along RA (x) 
and Dec (y) are also listed in Table HI The proper motion vector is listed in Table El as are 
astrometric results for the other PNNi, including catalog statistics. 

To determine if there might be unmodeled - but possibly correctable - systematic effects 
at the 1 mas level, we plotted the DeHt 5 reference frame x and y residuals against a number 
of spacecraft, instrumental, and astronomical parameters. These included x,y position within 
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the pickle-shaped FGS field of regard; radial distance from the center of the FGS field of 
regard; reference star V magnitude and B-V color; and epoch of observation. We saw no 
obvious trends, other than an expected increase in positional uncertainty with reference star 
magnitude. 

5.4. The Absolute Parallax of the DeHt 5 Central Star 

For this object at high ecliptic latitude (note the large parallax factors in both RA and 
Dec in Table [2]) we can solve for the separate x and y components of the parallax. These 
were vr^ = 2.79 ± 0.17 mas and iTy = 3.19 ± 0.35 mas. We obtain for the DeHt 5 PNN a 
final absolute parallax T^abs = 2.86 ± 0.16 mas. Our result agrees within the errors with 
the previous ground-based parallax measurement of the DeHt 5 PNN (Harris et al. 2007), 
TTafes = 3.34 ± 0.56 mas. Parallaxes from HST and USNO and relative proper motion results 
from HST are collected in Table O Even though both proper motion determinations are 
relative, using different sets of reference stars, the proper motion vector agrees with that 
determined by USNO (Harris et al. 2007). For the remainder of this paper we adopt as 
the absolute parallax of the DeHt 5 PNN, Habs = 2.90 ± 0.15 mas, the weighted average of 
these two independent parallax determinations. The degree of independence is quantified in 
Table O Of the ten reference stars used in the HST study, only three were in the suite of 
fifteen (on average much fainter) reference stars used in the USNO study. If both studies 
used exactly the same set of reference stars, one would expect that some component of the 
uncertainties would be correlated. As indicated in Table [5l this is not the case for any of the 
PNNi investigated in our study. 

5.5. Modeling Notes on the Other PNNi 

Abell31 - This field provided six reference stars. The reference star average data are 
listed in Table [51 We again used the eight parameter model (Equations 4 and 5). Because of 
the low ecliptic latitude, most of the parallax signature is along RA. Hence, we constrained 
= T^y Two gyro guiding scheduling constraints and the afore-mentioned science-side 
problems yielded a total study duration of four years. The HST parallax, iTabs = 1-51 ± 0.26 
mas agreed well with the USNO value from Harris et al. (2007), TCabs = 1-76 ± 0.33. For 
the remainder of this paper we adopt as the absolute parallax of the Abell31 PNN, iTabs = 
1.61 ± 0.21 mas, the weighted average of these two independent parallax determinations. 
HST , USNO, and final HST + USNO weighted average parallaxes are given in Table O 
This distance, d= 621^70 PC, rules out a physical association between DeHt 5 and the red 
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companion detected by the HST WFPC-2 camera discussed in ICiarduUo et al.l (119991 ) 



NGC 7293 - This field provided only three useful reference stars. The reference star 
average data are listed in Table O Because of the paucity of reference stars, the astrometric 
model for this field used only four parameters, discarding the radial terms and constraining 
d=-b, e=a in Equations 4 and 5. We also constrained Hx = iiy. Model selection was dictated 
by the loss of access to 1-2 reference stars for many of the observation sets, primarily due to 
two-gyro guiding constraints on allowed spacecraft roll. One of the reference stars (#18 in 
the original GO-10432 proposal) was removed from our modeling. For several observation 
sets each visit to that reference star locked on a different component. From these systematic 
residuals after initial modeling we inferred that it is either a binary or an optical double with 
component separation of ~ 11 mas. Our NGC 7293 PNN parallax is Hats = 4.67 ±0.16 mas. 



van LeeuwenI (120071 ) asserts th at the intrinsic width of the main sequence for FGK 
stars is 0.4 mag (l-cx). Past results (IBenedict et al.l 120071 ) indicate that for fields with 5 or 
more reference stars, cosmic dispersion in reference star absolute magnitude has no apparent 
consequence. In fact, Cepheid astrophysics argues that our parallax errors are overstated. 
To obtain a unity reduced ix^ I DOF where DOF = degrees of freedom) for our linear 
My - logP Period-Luminosity relation, we must significantly reduce our magnitude errors. 
For this one target whose parallax is dependent on only three reference stars we explored 
the effects of cosmic dispersion on reference star input absolute parallaxes. Worst-case (l-a 
increase or decrease in the My of all reference stars) the final absolute parallax for NGC 7293 
could range ±0.29 mas. Hence we add in quadrature that error due to cosmic disper s ion to 
our astrometry-only error for a final result of i^ahs = 4.67 ± 0.33 mas. [Harris et al.l (120071 ) 
obtained Habs = 4.56 ± 0.49 mas. A weighted average yields iiabs = 4.66 ± 0.27 mas. Our 
measured proper motion vector does not agree within the errors with the USNO value. This 
is not unexpected, given that our proper motion is measured against so few reference stars. 



NGC 6853 - As noted in [Benedict et al.l ( 120031 ) our original FGS 3-only data did not 
adequately cover the epochs of maximum parallax factor, resulting in a relatively (for HST ) 
poorly determined parallax, TCats = 2.10 ± 0.48 mas. The addition of two new epochs of ob- 
servation at maximum parallax factor with FGS Ir significantly increased the duration of the 
study (now over nine years) and improved the precision and slightly increased the parallax, 
now TTabs = 2.22 ± 0.19 mas. The astrometric model for this field used only six parameters, 
discarding the radial terms in Equations 4 and 5, a choice dictated by an insignificant decline 
in reduced when increasing from 6 to 8 astrometric coefficients (Equations 4 and 5). We 
also constrained vr^. = iTy. For this object there are two USNO results (Harris et al. 1997, 
Harris et al. 2007) which yield the weighted USNO average Hats = 3.17 ± 0.32 mas and the 
final HST + USNO weighted average iTabs = 2.47 ± 0.16 mas listed in Table [51 We note the 
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significant difference in relative proper motion between HST and USNO. 



6. PNNi Absolute Magnitudes, Radii, and Masses 

Again, we use DeHt 5 as an example of the steps required to obtain absolute magnitudes 
and radii for these PNNi. The final results for the other three PNNi are summarized in 
Table [6] and in individual notes (Section 16.41) below. 



6.1. Absolute Magnitudes and the Lutz-Kelker-Hanson Bias 



When using a trigonometric parallax to estimate the absolute magni tude of a star, a 
correction sho uld b e made for the Lutz-K elker bias (ILutz fc Kelken Il973l ) as modified by 
HansonI (119791 ). See [Benedict et al.1 (120071 ). section 5, for a more detailed rationale for the 
application of this correction to single stars. Because of the galactic latitude and distance of 
DeHt 5, and the scale height of the stellar population of which it is a member, we calculate 
Lutz-Kelker-Hanson (LKH) bias twice, assuming first a spheroidal then a disk distribution. 
The LKH bias is proportional to (cTtt/tt)^. Presuming that the PNN belongs to the same class 
of object as 6 Cep (young, evolved Main Sequence stars in a core helium burning phase), we 
scale the LKH correction determined for S Cep in Benedict et al. (2002b) and obtain LKH 
= -0.02 magnitude. Presuming that the PNN belongs to the same class of object as RR Lyr 
(older, evolved Main Sequence stars on the horizontal branch), we scale the LKH correction 
determined for RR Lyr in Benedict et al. (2002a) and obtain LKH = -0.03 magnitude. 
Our final LKH bias corrections are an average of the biases from the two adopted prior 
distributions. The corrections differ by at most 0.02 magnitude. See Benedict et al. (2007), 
section 5, for a more detailed rationale, justifying the use of this correction to single stars. 



6.2. The Absolute Magnitude of the PNN of DeHt 5 

Adopting for the DeHt 5 PNN V= 15.47± 0.03 and the weighted average absolute 
parallax, TTabs = 2.90 ± 0.15 mas from Table [5], we determine a distance modulus, m-M = 
7.69±0.12. To obtain a final absolute magnitude, we must correct for interstellar extinction. 
There are a variety of techniques used to estimate the extinction towards (and internal to) 
planetary nebulae. Perhaps the most common is the assumption that in nebular conditions, 
the ratio of Ha/H/3 = 2.86, and any deviation from this value is assigned to extinction. 
Alternatively, the observed radio flux can lead to a prediction of the flux in H/3 assuming an 
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optically thin nebula with a temperature of 10 K. Either technique leads t o a value for a 



logar ithmic extinction at H/? usually denoted by "c" , where E{B — V) = 0.83c (IMilne fc AUer 



19751 ). Other techniques include assuming a single value for the color for the PN central star, 
using interstellar absorption features seen in the optical or ultraviolet spectra, or simply using 
field stars along similar lines of sight. 

For DeHt 5 tw o estimates of "c" (Pottasch 1996, Phillips 2005) hst E{B - V) = 0.0. 



Harris et all (120071 ) derive E{B - V) = 0.18, assuming {B - V)o = -0.38, yielding (R = 
Ay/E(B-V) = 3.1) = 0.56. Estimating with field stars, from Table E] (Section USD we 
derive a per-star, per-unit 100 pc distance absorption, (Ay) /lOOpc. The average of the three 
stars nearest the central target (see Figured]), ref-2, -4, and -6, is (Av)/100pc = 0.09±0.01. 
With this per-unit 100 pc (Ay) and the distance to the DeHt 5 central star, d — 345^^^ pc, 
we obtain a total absorption for the PNN, Ay = 0.32 ±0.03. We also estimat e Ay us i ng th e 



measured temperature and a grid of synthetic spectra of hot, compact stars iRauchI (120031 ). 
We calculate an intrinsic {B — V)o = —0.36. This yields Ay = 0.42 ± 0.07, where the error 
is obtained through a 50000 trial Monte-Carlo process. Given the scatter in Ay from the 
various determinations, we choose to average the determination from synthetic spectra and 
the reference star values, yielding Ay = 0.37 ± 0.07, deeming the 'c' determination fiawed 
for this faint an object. Including the LKH correction we obtain an absorption-corrected 
magnitude, Vq = 15.03. The distance modulus and Vq provide an absolute magnitude 
My = 7.39l°:it 



6.3. A Radius for the PNN of DeHt 5 

We employ two methods to derive a radius, both differential in nature. The first method 
employs the Stefan-Boltzmann relation, the second filter-averaged Eddington fiuxes. 

To estimate a radius, R^,, for this star using the Stefan-Boltzmann relation we require a 
distance, an absolute magnitude, an effective temperature, T^ljf, and a bolometric correction 
(B.C.) . These quantities then yield a radius via differential comparison with the Sun. Our 
parallax provides a distance, d = 345^^7 pc and an absolute magnitude. My = 7.39 ± 0.14. 



Napiwotzkil (119991 ) has estimated T*jj = 76, 500K± 5800 K from model atmosphere fits 
to the Balmer H5 and He absorption lines. We calculated bolometric corrections from our 
synthetic photometry, by comparing the integrated surface fiux crTjj with the surface flux 
through the V band fllter including the filter constant (or in other words computing the offset 
between Mboi and My of a star with arbitrary radius). We set the zeropoint by adopting 
the following values for the Sun: = 4.75 and My = —26.74. For DeHt 5 we calculate 

B.C. = - 5.84 ± 0.2, where the error is dominated by the uncertainty in Tiff- 
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We obtain a PNN bolometric luminosity Mboi = My + B.C. = +1.55 ±0.24. follows 
from the expression 

- M,;, = 10 logiT:^^/ T®^) + 5 log{Rj R^) (6) 

where we assume for the Sun Tf^^ = 5800i^. We find -R* = 0.026 ± 0.007/?©. The sources 
of error for this radius are in the absolute magnitude (i.e., the parallax), the bolometric 
correction, and the ^e*//- 



m 



A second way to ob tain R^ involves the V-band aver age Eddington flux, Hy, discussed 



Bergeron et al.l (1l995l ) and iHolberg &: BergeronI (l2006l ) . The latter carried out a careful 



photometric calibration of DA white dwarfs with pure hydrogen atmospheres. However, 
radiative levitation in the hottest white dwarfs causes metals to be present in the atmospheres 
on a level roughly equivalent to solar abundances (Barstow et al. 2003). This has some 
impact on the SED causing redistribution of flux from the UV to longer wavelengths. To 
take this into account we computed synthetic photometry in the Johnson system from NLTE 
model spectra calculated with solar abundances of important elements up to the iron group 



(IRauchl l2003l ). The photometry is linked to the Vega system as outlined by Holberg & 
Bergeron (2006). The overall effect of the metals is that B — V colours for hot PNNi are 
bluer by ~0.03 mag, and the flux level in the V band increases by typically 0.2 mag compared 
to pure hydrogen atmospheres. 

Our synthetic photometry provides Hy as a function of temperature for solar-met allicity 
WD of various temperatures calibrated against Vega. We obtain Hy for TJ^-j = 76, 500 K. 
With H^''", we can derive i?* from 

Rl = (ffj:-57i7*)io-o-4(^v-<^''") (7) 

We obtain for DeHt5 with T*^^ = 76,5007^, = 7.39 from our parallax, and mJ^'^" = 
+0.026 an i?* = 0.025 ± 0.002/?©, where the radius error is obtained through a 50000 trial 
Monte-Carlo process. Given that the approach relying directly on the B.C. and the approach 
utilizing Hy yield R^, values that agree, we will use for the remainder of this paper the lower 
error, R^, = 0.025 ± 0.002/?©. For the other three objects the radius from Stefan-Boltzmann 
will be calculated as a confirmation only. The higher error from Stefan-Boltzmann is due 
primarily to the significant contribution to the error budget from the bolometric correction 
uncertainty. The error on this radius cannot be further reduced by a weighted average of 
the results from the two approaches, because their errors are highly correlated, both having 
contributions from the uncertainties in T*jrj and My. 
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6.4. Notes on the Radii of Other PNNi 



All quantities required to estimate the radii for our other PNNi, as we did for DeHt 5, 
are gathered in Table [6l where we also list our derived radii. 



AbellSl - iKalerl (119831) estimates V = 0.0 ±0.4. We note that the line of sight total 
extinction estimated by ISchlegel et al.l (119981 ) is E{B ~ V) = 0.064. In this case we derive 
an average absorption per 100 pc for the four astro metric referen ce stars nearest in angular 
distance to the PNN of Ay = 0.05 ± 0.04. From iRauchI (120031 ) we calculate an intrinsic 
(B - 1/)o = -0.36. This yields k*y = 0.24 ± 0.07. The average of these methods yields 
Ay = 0.10 ± 0.07 which we adopt. The distance modulus and Vq provide an absolute 
magnitude My = 6.31^q;26, where we have included an LKH correction (-0.14 ± 0.03) and 
its uncertainty and the small uncertainty in in quadrature. 



Napiwotzkil (119991 ) has estimated T*^^ = 84, 700K± 4,700 K. From our synthetic pho- 
tometry we calculate a B.C. = - 6.29 ± 0.2. As in Section 1^751 above, we compare bolometric 
luminosities with the Sun and find = 0.041 ± 0.017i?f7,. We obtain H?, for T* 



K, again from our synthetic photometry, which yields -R* = 0.039 ± 0.006i?, 



■eff 



84, 700 



■0- 



NGC 7293 - has a reliably low value for the extinction. iKalerl (1l983l ) fin ds a value of 



Bohlin et al. 



c = 04. Milne & AUer (1975) derive E{B - V) = 0.01 from B.0 and radio . 
(Il982h using lUE observations derives E(B-V)=0.012+/-0.03. Iflarris et all J2007h quotes 
E(B-V)=0.03 from BVI photometry of the central star (assuming [B — V)o = —0.38), while 
Pottasch (1996) u s es the sam e technique and gets E{B — V) = 0.0. The extinction maps of 
Burstein fc HeilesI (119821 ) and lSchlegel et al.l (119981 ) would estimate a red dening of E(-B — V) 
~ 0.0 3 mag, consistent with the Nal measurements for the central star (IMauron &: Kendall 



20041 ). If we translate the values for c to E{B — V), and incorporate all of the other 
values we get a mean of {E{B — V)) = 0.027 ± 0.022. Consistent with this value, we 
derive an average absorption per 100 pc from the three astrometric reference stars nearest 
in angular distance to the PNN. With a per-unit 100 pc (Av)=0.02±0.02 and the measured 
distance to the NGC 7293 central star, d = 216±13 pc, we obtain a total absorption for the 
PNN, A;^ = 0.04 ± 0.03. Again we adopt an average of the estimates, A;^ = 0.09 ± 0.04, 
yielding My = 6.74 ± 0.13, where we have included an LKH correction (-0.03 ± 0.01) and 
its uncertainty and the small uncertainty in Ay in quadrature. 



Napiwotzkil (119991 ) has estimated T* 



// 



103, 600K± 5,500 K, which yields B.C. = -6.77. 
Comparing bolometric luminosities with the Sun we obtain R^, = O.O28±O.OO7i?0. We obtain 
103, 600 K from our synthetic photometry, which yields i?* = O.O28±O.OO3i?0. 



for T:^^ 



NG C 6853 - To es t imate A,/ iBarkeil (Il984[ ) found c = 0.17 and lKaler et all (119761 ) find 
0.02. iBohhn et all Jl982h estimate E(B-V)=0.06+/-0.03 from lUE observations, while 
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Pottasch (1996) gets E( B - V) = 0.10 andlHarris et al.l (120071 ) find E(B - V) = 0.07 from 



central star pliotometry. ICiardullo et al.l (119991 ) obtain c=0.11. The mean with error is then 
{E{B - V)) = 0.08 ± 0.06, or (Ay)=0.26 ± 0.17. Given this uncertain (Ay), we derive an 
average absorption per 100 pc f or the three star s near est the central target, ref-4, -5, and -8. 
A finder chart can be found in [Benedict et al.l (120031 ). With a per-unit 100 pc (Ay)=0.07 
± 0.03 from these three stars and the measured distance to the NGC6853 central star, d 
= 405+25 pc, we obtain a total absorption for the PNN, A;^ = 0.28 ± 0.06. We adopt an 
average Ay = 0.30 ± 0.06 and obtain My = 5.62 ± 0.16, where we have included an LKH 
correction (-0.03 ± 0.01) and its uncertainty and the 0.06 magnitude uncertainty in Ay in 
quadrature. 



Napiwotzkil (119991 ) has estimated T*^^ = 108, 60QK± 6800 K for which we obtain B.C. 



= -6.91 ± 0.2, where the error is dominated by the uncertainty in the temperature and 
the behavior of the B.C. at these high temperatures. Comparing bolometric luminosities 
with the Sun we find i?* = 0.046 ± 0.0127?©. Comparing Hy with Vega we obtain Hy for 
T*fj = 108, 600 K hsted in Table El which then yields = 0.045 ± 0.004i?Q. Given that the 
two approaches yield R^, values that agree within their errors, we adopt i?* = O.O45±O.OO4i?0. 



6.5. Astrophysical Consequences 



6.5.1. Radii: PNNi vs. WDs 



Our four parallaxes, along with measured temperatures and apparent luminosities have 
resulted in four newly estimated radii for PNNi that, according to theory, should eventually 
descend to a WD cooling track. Previous investigations have yielded precise temperatures 
and radii (and masses) for five WD in visual, spectroscopic, and eclipsing binaries. These 
results (Sirius B, Holber g et al. 1998; Procyon B , Girard et al. 2000; 40 Eri B, Shipman et 



al. 1 997) are collected in jProvencal et al.l (120021 ). To these we add Feige 24 (IBenedict et al. 
20001). which has (T ^/tt = 2.7%, hence, a reasonably well-determined radius, two PNNi from 



Harris et all (120071 ). Sh 2-216 and HDW 4, b oth with ct^/tt < 8%, and the eclipsing Hyades 
WD binary, V471 Tau (lO'Brien et al.l 1200 ll ) . We will now compare PNNi and WD, using 
the quantities gathered in Table [71 We can also test the accuracy of the parallax and 
the many corrections leading to the bolometric magnitudes required by Equation 6, relating 
temperature, bolometric magnitude and radius. This aggregate of data probes the transition 
from PNNi to WD. 



Let us for the moment assume that all PNNi have the same mass and radius. In such 
a Universe the only PNNi and/or WD variable is age, hence, temperature, and there should 
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be a relationship between absolute bolometric magnitude and temperature, through Stefan- 
Boltzmann, L=4n Rl^aT^. Figure [7] (basically an HR diagram) indicates that temperature 
and absolute bolometric magnitude are correlated, in that hotter PNNi and WD have brighter 
absolute magnitudes. This is not surprising, given that the measured or estimated mass of 
most of these objects lie in the range 0.44 < J^p) < 0.60. For the WD, bolometric magnitudes 
are derived using B.C. from Flower (1996) and lBergeron et al.l ( 119951 ). The most discrepant 
object is Sirius B, which has the greatest mass, A4* = l.O2A^0, hence, the smallest WD 
radius. If Equation 6 perfectly describes both PNNi and WD, then the residuals, AMhoi, 
in Figure [7] to our simple linear relationship between Mtoi and log T should correlate with 
the log of the radius. Figure [H wherein AM^oi is plotted against log radius, shows such a 
correlation. The residuals in Figure[8]have an RMS dispersion of 0.12 magnitudes, indicating 
the accuracy of the LKH and extinction corrections. The two most discrepant objects are 
Feige24 and Sirius B with residuals of order 0.2 magnitudes. 



6.5.2. PNNi Masses and Gravities 

Various investigators have modeled the evolution of a star as it passes through the 
red giant phase, ejects significant mass, and becomes a white dwarf. The PN phase lies 
between the giant and WD stages in stellar evolution. To estimate PNNi mass we compare 
in Figure [3 the positions of our PNNi in an H-R diagram (My - log Tg//) with predicted 
evolutionary tracks of post-AGB sta rs from Schonberner & Blocker (1996). Also plotted 
are tracks of lower mass stars from iDriebe et al.l (119991 ). Absolute V-band magnitudes. 
My, for tracks were c alculated for solar metallicity, using SEDs from Ranch (2003). PNNi 
temperatures are from iNapiwotzkil (119991 ). We find that the four PNNi clump about a mass 
M = 0.57M 0, in agreement with the peak of the WD mass distribution, 0.60 A4q, found 
by Liebert, Bergeron & Holberg (2005). Interpolated individual masses are found in Table [61 
We also estimate masses for PNNi HDW4 and Sh2-216. These are listed in Table [3 We 
note that the mass of 40 Eri B from these tracks is ~ O.9A^0, differing substantially from 



that listed in jProvencal et al.l (Il998l ) 



A mass and radius uniquely determine a surface gravity, g, through 

g = MG/R^ 



where G is the gravitational constant. Our DeHt5 radius, -R* = 0.025 ± O.OO2_R0 and the 
mass inferred from the evolutionary tracks in Figure [TOl (where we plot radius vs. Tg//), 
}A = 0.57 ± O.O2A^0, yield log g = 7.41±0.08. The uncertainty in our log g is primarily 
due to the radius uncertainty. For Abell31 our mass estimate, Ai = 0.53 ± 0.037^0 with 
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our radius determines a gravity, log g = 6.96±0.14. For NGC7293 using the radius from 
the Hy approach and a mass estimated from Figure [9|, A4 = 0.60 ± O.O2A^0 determines a 
gravity, log g = 7.34±0.06. For NGC6853 our mass estimate, M = 0.57 ±0.01A^,t^ w i th our 
radius determines a gravity, log g = 6.89±0.08, in agreement with the iNapiwotzkil (jl999l ) 
line profile fitting value, log g = 6.7ib0.2 . Calculated value s of log g for our four PNNi are 
listed in Table [6] and compared with the INapiwotzkil (119991 ) line profile fitting values. Note 
that Figure [10] shows 40 Eri B to have a mass consistent with past estimates. In Figure [11] 
we compare our WD sample with the PNNi listed in Table [7] on a mass-radius diagram. 
Most of the PNNi and the hottest WD , Feig e 24, lie above the zero-temperature mass-radius 
relationship from lHamada fc Salpeterl (Il96ll ). These radii confirm that PNNi are larger than 
cooler WD. 



6.5.3. Limits on Stellar Companions 



With the parallaxes from Table [S], absolute magnitudes and mass estimates from Table IS] 
and the separation and Am detection limits from Section 13.21 we can now estimate the 
spectral type, separation in AU and periods for companions at the limit of detectability. 
Companions with spectral types later than listed in Table [8] would not be detected by the 
FGS. For example a companion two magnitudes fainter than the DeHtS PNN (an MIV 
star) is at the limit of detectability for a separation of 15 mas, which for the parallax of 
DeHt 5 equates to 5.2 AU. From P^{Mi + M2) = we derive a period P=lly. For another 



application of Table [8] we consider NGC7293. ISu et al.l (120071 ) find evidence for IR-emitting 
dust 35-150 AU from the PNN of NGC7293. A binary companion could provide an engine 
to sculpt a dust disk. Given that 35 AU at the distance of NGC7293 is 0'.'16 and that we 
would detect a companion with Am < 3.5, that companion would have to ha ve a spectral 



type later thaii MIV. If much asymmetrical PN structure is due to binarity (ISokerl 12006 



De Marcoll2009l ). then we would have the highest probability of detecting the companion to 
the PNN of NGC6853, the most asymmetric of the PN we observed. Our null detection 
suggests that such a companion is likely to have P<7y and M. < 0.7J^q. 

Finally, in principle we can probe for very short period companions using the photometry 
in Figure 2. Companions with small separations could produce a single-peaked orbital light 
curve through heating of the companion star by the PNN e.g., the reflection effect. As the 
companion star orbits the PNN, its heated face is alternately more or less visible, increasing 
and decreasing the observed flux from the PNN once per orbit. For example, the Feige 24 
system (WD+M2 V) has a period of 4.23*^ and f rom FGS pho t ometr y (Benedict et al. 2000) 
evidences a photometric variation of 25 mmag. iKawka et al.l (120081 ) estimate an inclination 
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i = 77°. Assuming a similar inclination, we find from binary light curve modeling (cf. 
Harrison et al. 2009) that even the shortest period companions would produce less V band 
variation from reflection effects than we see in Figure 2. PNNi are too bright, washing out 
any variations due to reflection from companions. 



6.5.4- Comparison with other Distance Estimates 



Napiwotzki (2001) has compiled PNe distances from a number of methods. Figure 
plots our distances from Table [6] against three other methods; distances obtained via non- 
LTE PNNi atmosphere analysis (Napiwotzki 2001), from an H/? -diameter relation devised by 
Shklo vski (1956), and from an interstellar Na D line analysis (INapiwotzki fc Schoenberner 
19951 ). Shklovskii distances are only available for only three of our targets, and the Na D 
method has been applied only to two o f our targets. The recent recalibration of the Shklovski 
distances by IStanghellini et al.l (120081 ) left those three distances basically unchanged. The 
dashed line represents perfect agreement. The Shklovski approach has a tendency to under- 
estimate the distances while the spectroscopic distances are a bit on the high side. We now 
provide a more detailed assessment of distances obtained via non-LTE PNNi atmosphere 
analysis. 



6.5.5. Reassessing the Spectroscopic Distance Scale for PNNi 

Napiwotzki (1999) determined PNNi fundamental stellar parameters, temperature and 
surface gravity, from a fit of the hydrogen Balmer lines with profiles computed from NLTE 
model atmospheres. This technique is well established and tested for the analysis of hot white 
dwarfs (e.g. Finley et al. 1997) using LTE atmospheres. However, when Napiwotzki (1999) 
applied this method to the even hotter central stars of old PNe, it became clear that for many 
stars no consistent fit of all Balmer lines could be achieved. A strong temperature trend was 
present, with the fit of higher members of the Balmer series yielding higher temperatures. 
This became known as the Balmer line problem (Napiwotzki & Ranch 1994). 



Napiwotzkil (119931 . Il999l ) presented arguments that the temperature derived from the 
highest Balmer lines B.5 and He are close to the real temperatures of the PNNi. However, 
a physical explanation of the Balmer line problem remained elusive for some time. Models 
used for the Napiwotzki (1999) investigation were calculated in full NLTE but included only 
the two most abundant elements hydrogen and helium. Tests carried out prior to the start 
of this project appeared to show that the impact of line blanketing of heavier elements on 
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the temperature structure of the atmospheres had only minor impact on the hydrogen hne 
profiles (see discussion in Werner 1996). However, Werner (1996) could show that strong 
cooling by the resonance lines of carbon, nitrogen, and oxygen can reproduce the observed 
Balmer lines in the hot sdO BD+28°4211 and the PNN of Sh 2-216, if detailed treatment of 
the Stark broadening of these lines is included in the calculations. This treatment is a very 
computer time intensive. For this reason it was not included in previous calculations. The 
Werner (1996) results provided an astrophysical explanation and essentially validated the 
recipe of adopting the temperature fitted for H5 and He put forward in Napiwotzki (1993). 
This method was then adapted by Napiwotzki (1999). 

Fitted surface gravities appeared to be unaffected by the Balmer line problem. All 
Balmer lines could be fitted with a single value of g using the Napiwotzki (1999) H and He 
models. Also, the Werner (1996) calculations did not indicate offsets in gravity. Napiwotzki 
(2001) compared his spectroscopic distance estimates with the results of other distance est i- 



mates including the best parallax measurements available at that time (IHarris et al.l 119971 ). 
The comparison showed spectroscopic distances larger than trigonometric distances by 55% 
- at face value. However, as pointed out by Napiwotzki (2001), one has to take into account 
that the often significant relative errors of trigonometric parallaxes introduce sample biases. 
Lutz-Kelker biases are one way to estimate the size of the effect. Na piwotzki (200 1 ) per- 



formed a Monte-Carlo simulation trying to model the properties of the [Harris et al.l (119971 ) 
sample as closely as possible given the - by necessity - not well defined selection criteria. The 
result was an estimated bias of the trigonometric distances, now too small by 32 ± 0.25%. 
The conclusion at that time was that both distance scales are marginally consistent, but 
large uncertainties remained. 

Improved accur acy of recent t rigon ometric parallax measurements have changed the sit- 



uation dramatically. [Harris et al.l (120071 ) achieved accuracies ranging from 0.3 mas to 0.6 mas. 



The measurements presented in our inv estigation yield even better accuracies with cTtt being 



0.21 mas and below. iHarris et al.l (120071 ) estimates a bias of 5% for their sample. A straight- 
forward reading from Table A. 2 in Napiwotzki (2001) gives an estimate of 7% for a sample 
with 0.2 mas accuracy. Both estimates confirm that remaining systematic errors are now 
much smaller, in line with the small Lutz-Kelker biases given in our Table Ei Thus it is a 
good time to re-assess the spectroscopic distance scale. 

Spectroscopic distances of the four program stars are 38% larger than the trigonometric 
distances (unweighted average). This offset is smaller than that found by Napiwotzki (2001), 
but due to the smaller errors and biases it is now highly significant. This translates into 
an average log (7 offset log ^fgpec — log (7,^ = —0.41. In Table [6] we find DeHt5 anomalous, 
the difference between gravity derived from our radius and from the analysis of stellar at- 
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mospheres larger than for the other three PNNi. DeHt 5 is of special interest and will be 
discussed below. Excluding DeHt 5 from the average we compute an average distance offset 
of 40% (spectroscopic distances recomputed using the improved photometric observations, 
reddening estimates and synthetic photometry) and \ogg offset of 0.30, not much larger than 
typical gravity errors given in Napiwotzki (1999) which are 0.2-0.3 dex. 



6.5.6. The case of DeHt 5 



The temperature and gravity derived by Napiwotzki (1999) place the central star of 
DeHt 5 in a region of the temperature gravity diagram inconsistent with a post-AGB origin. 
The parameters of this central star were better matched to those of a star which lost its 
envelope at the end of the first red giant branch (RGB) and is now evolving into a low 
mass He-core white dwarf. iBarstow et al.l (120011 ) performed an analysis of optical and UV 
spectra using model atmospheres including the effect of model line blanketing. The derived 
temperature is 57,400 K - lower than the Napiwotzki (1999) result - partly due to the 
inclusion of metal line blanketin g and partly due to a different fit algorithm and philosophy. 
The gravity log (7 = 7.0 given by iBarstow et al.l (1200 ll ) is essentially an upper limit, because 
it was the lowest log (7 available in the model atmosphere grid. The gravity resulting from 
a fit with pure hydrogen models is in good agreement with the Napiwotzki (1999) result: 
log^f = 6.75 vs. 6.65. 

None of the results can be recon ciled with the trigonometric results, which translates 
into log (7 = 7.4. IBarstow et al.l (12001! ) also determined metal abundances from the analysis 
of an HST-STIS UV spectrum. The resulting abundances are higher than those of the well- 
known hot "template" white dwarf G191-B2B, which appears to have typical abundances 
for that parameter range (Barstow et al. 2003). One could be tempted to speculate that 
unusual metal abundances could explain the unusual large differences between trigonometric 
and spectral analysis results. More detailed abundance analyses of more PNNi would be 
needed to decid e this question . In an y case one conclusion is that even the metal blanketed 
atmospheres of IBarstow et al.l (120011 ) are not capable of producing results consistent with 
the trigonometric parallaxes. 

The stellar parameters implied by the trigonometric parallax {R = O.O25i?0 or logg = 
7.4) places DeHt 5 at a location expected for a run-of-the-mill pre-white dwarf. The implied 
mass M = 0.57Mq sits spot on the main peak of the white dwarf mass distribution (M = 
0.572; Liebert et al. 2005). The implication is that DeHt 5 is a rather normal C/0 white 
dwarf resulting from post-AGB evolution. Somewhat problematic is the high implied post- 
AGB age of > 4 X 10^ yr (read from the O.6O5M0 track). The observational sample implies 
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that PNe disperse into the intersteller medium after about 50 — 100 x 10^ yr (e.g. Table 1 in 
Napiwotzki 2001). A post-AGB age as high as imphed for DeHtS would be highly unusual. 
Parker et al. (j2006 ) speculate that DeHtS is not a real PNe, but a chance association of an 



interstellar cloud with a hot white dwarf. A more detailed investigation of the nebula will 
clarify this issue. 



7. Summary 

HST FGS photometry indicates that none of these PNNi shows photometric variation 
larger than 5 mmag. From FGS interferometric fringe morphology we establish companion 
limits mid-KV to early MV. FGS interferometric fringe tracking astrometry yields abso- 
lute trigonometric parallaxes for the PNNi of DeHtS, Abell31, NGC7293, and N GC68S3. 



Weighted averages with previous ground-based determinations (IHarris et al.l 120071 ) provide 
parallaxes with errors at or below 0.2 mas, or (a^/vr) = S%. Our results confirm that statis- 
tical distances methods of the Shklovski type underestimate the distances of old planetary 
nebulae. On the other hand, the improved accuracy of our trigonometric parallaxes now 
show that previous spectroscopic distances significantly overestimated the distances. We see 
a consistent trend in the spectroscopic distance scale overestimating the true distance by 
40% (corresponding to an understimate of log (7 by 0.3 dex). Results from Napiwotzki (1999) 
and similar studies should be corrected accordingly. We use these parallaxes and estimates 
of interstellar extinction from our own spectrophotometry and other investigations to de- 
rive PNNi absolute magnitudes. With these we derive radii, either comparing with the Sun 
through bolometric corrections, or with Vega via the V-band average flux, Hy- These four 
PNNi along with two others with well-determined distances and five WD satisfy theoretical 
linear correlations between absolute bolometric magnitude, log temperature, and log radius. 
Estimating from post-AGB evolutionary models, we find PNNi masses that agree with those 
typically found for white dwarf stars. The PNNi and the hottest WD clearly fall above a 
WD mass-radius relationship established by nearby, cool WD. 
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Table 1. PNNi Positions and Aliases 



FN 


RA (2000) 


Dec 




Aliases 




NGC6853 


19 59 36.34 


+22 43 16.1 


Dumbbell 


M 27 


WD 1957+225 


NGC 7293 


22 29 38.55 


-20 50 13.6 


Helix 


FN G036.1-57.1 


WD 2226-210 


DeHt5 


22 19 33.71 


+70 56 03.3 




FN Gill. 0+11.6 


WD 2218+706 


Abell 31 


08 54 13.16 


+08 53 53.0 


FN A66 31 


FK 219+31 


FN CSI+09-08515 



Table 2. DeHt 5 Log of Observations and Reference Star Availability (x = Observed) 



Set 


mJD 


p a 
^ a 




DeHt5 


2" 


3 


4 


5 


6 


7 


8 


9 


10 


11 


1 


53390.0556 


-0.5020 


-0.8415 


X 


X 


X 


X 


X 




X 


X 


X 


X 


X 


2 


53392.0128 


-0.4752 


-0.8578 


X 


X 


X 


X 


X 




X 


X 


X 


X 


X 


3 


53587.7651 


0.3155 


0.9660 


X 


X 


X 


X 


X 




X 


X 


X 


X 


X 


4 


53587.8598 


0.3147 


0.9663 


X 


X 


X 


X 


X 




X 


X 


X 


X 


X 


5 


53623.0390 


-0.2369 


0.9536 


X 


X 


X 




X 


X 












6 


53671.1814 


-0.8217 


0.4197 


X 


X 


X 


X 




X 












7 


53764.2406 


-0.3764 


-0.9099 


X 


X 


X 


X 


X 




X 


X 


X 


X 




8 


53772.1216 


-0.2563 


-0.9501 


X 


X 


X 




X 




X 


X 


X 


X 


X 


9 


53783.0153 


-0.0818 


-0.9759 


X 


X 






X 


X 












10 


53956.0037 


0.2697 


0.9764 


X 


X 


X 


X 


X 


X 












11 


53957.1156 


0.2521 


0.9805 


X 


X 


X 


X 


X 


X 












12 


54780.6123 


-0.8997 


0.1901 


X 


X 


X 


X 




X 












13 


54780.6788 


-0.9000 


0.1890 


X 


X 


X 


X 




X 












14 


54780.7454 


-0.9003 


0.1879 


X 


X 


X 


X 




X 












15 


54782.5436 


-0.9067 


0.1574 


X 


X 


X 


X 




X 













'^Farallax factor in Right Ascension 
''Farallax factor in Declination 

'^Reference star number 
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Table 3. DeHt 5 Astrometric Reference Star Photometry, Spectral Classifications, and 

Spcctrophotomctric Parallaxes 



ID 


V 


B-V 


V-K 


SpT 




Ay 


TTabs (mas) 


Ref-2 


15.56 


0.82 


2.29 


F5V 


3.34 


1.33 


0.66±0.15 


Ref-3 


14.90 


0.82 


2.2 


F4V 


3.11 


1.33 


0.77±0.19 


Ref-4 


14.31 


0.96 


2.54 


G2V 


4.56 


1.25 


1.98±0.46 


Ref-5 


11.92 


0.48 


1.24 


F4V 


3.11 


0.25 


1.84±0.45 


Ref-6 


13.55 


0.6 


1.63 


F4V 


3.11 


0.72 


1.07±0.26 


Ref-7 


15.02 


1.06 


2.76 


G2V 


4.56 


1.49 


1.60±0.37 


Ref-8 


14.90 


0.83 


2.72 


F2V 


2.84 


1.32 


0.85±0.17 


Ref-9 


13.50 


1.42 


3.75 


Kim 


0.6 


1.34 


0.49±0.12 


Ref-10 


14.94 


0.92 


2.52 


F7V 


3.72 


1.27 


1.05±0.24 


Ref-11 


14.86 


1.48 


3.88 


K2 III 


0.5 


1.34 


0.25±0.06 
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Table 4. DeHt 5 and Reference Star Astrometric Data 



ID 




77 ^ 








DeHt 5*^ 


-4.5797±0.0001 


-2.5814±0.0001 


-11.80±0.10 


-18.49±0.08 


2.86±0.16 


Ref- 


■2 


-3.5388±0.0002 


52.3476±0.0001 


-3.19±0.11 


-5.59±0.10 


0.70±0.05 


Ref- 


3 


-2.0182±0.0002 


-156.6239±0.0002 


-2.51±0.12 


-5.21±0.13 


0.90±0.03 


Ref- 


■4 


-92.2348±0.0003 


-92.3332±0.0002 


1.73±0.12 


1.77±0.16 


1.30±0.04 


Ref- 


-5 


53.3040dz0.0001 


161.1489±0.0001 


0.20±0.24 


-7.28±0.28 


1.50±0.05 


Ref- 


■6 


117.4994±0.0002 


-52.1638±0.0002 


2.15±0.10 


-4.14±0.10 


1.10±0.08 


Ref- 


■7 


-116.1855±0.0004 


-189.3931±0.0004 


6.41±0.55 


5.39±0.64 


0.80±0.15 


Ref- 


■8 


-132.3795±0.0003 


-311.4981±0.0004 


0.43±0.50 


0.10±0.47 


1.21±0.10 


Ref- 


■9 


-47.7228±0.0003 


-263.5350±0.0003 


-3.29±0.50 


-1.06±0.49 


0.34±0.16 


Ref- 


■10 


73.5099±0.0004 


89.0934±0.0003 


0.30±0.41 


1.17±0.51 


1.01±0.08 


Ref- 


11 


-8.2819±0.0006 


202.4660±0.0007 


-1.16±0.96 


1.07±1.02 


0.21±0.16 



(RA) and 77 (Dec) are relative positions in arcseconds. 
and Hy are relative motions in mas yr~^, where x and y are aligned with RA and Dec. 
■^Absolute parallax in mas 

^RA = 22'^19"'33.713" +70°56'03.28", J2000, epoch = mJD 53764.24692 
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Table 5. Reference Frame Statistics and PNNi Parallax and Proper Motion 



Parameter 



PNNi 





DcHt 5 


Abell31 


NGC 7293 


NGC 6853 


HST Study Duration (y) 


3.81 


3.99 


1.84 


9.10 


Observation Sets (#) 


15 


15 


11 


12 


Ref stars (#) 


10 


6 


3 


7 


Ref stars (V) 


14.29 


13.56 


13.11 


14.37 


Ref stars (B-V) 


0.94 


0.80 


0.69 


1.28 


< > (mas) 


0.3 


0.3 


0.9 


0.3 


< cr,, > (mas) 


0.3 


0.2 


0.9 


0.3 


HST Tiahs (mas) 


2.86±0.16 


1.51±0.26 


4.67±0.33 


2.22±0.19 


il/^ST Relative /i (mas y~^) 


21.93±0.12 


10.49±0.13 


38.99±0.24 


8.70±0.11 


in Position Angle (°) 


212.5±0.10 


227.1±0.2 


100.1±0.3 


67.9±0.11 



USNO Ref stars (# 
Ref stars (#) in common 
USNO TTabs (mas 
USNO Relative fi (mas y~^ 
in Position Angle (° 



15 
3 

3.34±0.56 

21.40±0.20 

214.6±0.5 



5 


1.76±0.33 

10.40±0.10 
226.5±0.6 



6 
1 

4.56±0.49 

33.0±0.1 

86.7±0.3 



28 
2 

3.17±0.32 

13.50±0.25 
60.8±1.0 



Weighted HST + USNO nabs 2.90±0.15 1.61±0.21 4.64±0.27 2.47±0.16 
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Table 6. PNNi Astropliysical Quantities 



Parameter 



PNNi 





T\ TTj_ f 

DeHt 5 


A 1 n o -1 

Abell 31 


NGC 7293 


NGC 6853 




15.47 


15.52 


13.53 


13.99 




-0.22 


-0.29 


-0.32 


-0.30 


d (pc)^ 


3451J? 


62i;?J 




405^1 


4* 


0.37 ±0.07 


0.10 ±0.07 


0.09 ±0.04 


0.30 ±0.06 


m-M 


7.69 ±0.12 


8.97 ±0.28 


6.67 ±0.13 


8.04 ±0.14 


LKH Bias 


-0.02 ± 0.02 


-0.14 ± 0.03 


-0.03 ± 0.01 


-0.03 ± 0.01 


Mv 


7.39 ± 0.14 


6.31 ±0.30 


6.74 ±0.13 


5.62 ± 0.16 




76, 500 ±5, 800 


84, 700 ± 4, 700 


103, 600 ±5, 500 


108, 600 ±6, 800 


B.C. 


-5.84 ±0.2 


-6.29 ±0.2 


-6.77 ±0.2 


-6.91 ±0.2 


M* 


±1.55 ±0.24 


±0.03 ±0.43 


-0.03 ±0.24 


-1.29 ±0.25 




1.12x10^ 


1.24x10^ 


1.62x10* 


1.72x10* 




0.025 ±0.002 


0.039 ±0.006 


0.028 ±0.003 


0.045 ± 0.004 




0.57 ±0.02 


0.53 ±0.03 


0.60 ±0.02 


0.57±0.01 


log g"" 


6.7 ± 0.2 


6.6 ± 0.3 


7.0 ± 0.2 


6.7 ± 0.2 


log 


7.41 ± 0.08 


6.99 ± 0.14 


7.34 ± 0.06 


6.89 ± 0.08 



''from [Harris et all (120071 ) 

'^from weighted average of HST and [Harris et al.l (120071 ) 




"ergs cm " s ^ A ^ str ^ 
■^from Fig. M 
Hrom g = MG/R^ 
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Table 7. Comparing PNNi and WD 



ID 


Vo 


(m-M)o 


My 


BC 


IV /T 

Mboi 


log T (K) 


log R*(©) 


IV /T / \ fJ" 

M(©)s 


DeHt 5^ 


15.15 


7.69 


7.39 


-5.95 


1.44±0.28 


4.884±0.04 


-1.61±0.03 


0.57±0.02 


AbellSl^ 


15.47 


8.97 


6.25 


-6.3 


-0.06 0.43 


4.928 0.03 


-1.40 0.06 


0.53 0.03 


NGC 7293^^ 


13.50 


6.67 


6.74 


-6.77 


-0.03 0.24 


5.015 0.03 


-1.56 0.06 


0.60 0.02 


NGC6853'^ 


13.75 


8.04 


5.62 


-7.09 


-1.47 0.25 


5.036 0.03 


-1.35 0.04 


0.57 0.01 


Sh 2-216^ 


12.38 


5.54 


6.82 


-6.30 


+0.52 0.23 


4.920 0.03 


-1.47 0.10 


0.55 0.03 


HDW 


16.11 


6.60 


9.45 


-4.47 


4.98 0.28 


4.674 0.03 


-1.87 0.12 


0.77 0.07 


Feige 24"= 


12.56 


4.17 


8.39 


-4.82 


3.57 0.13 


4.751 0.02 


-1.73 0.02 


0.57 0.02 


V471 Tau^ 


13.72 


3.39 


10.33 


-3.49 


6.84 0.03 


4.538 0.01 


-1.97 0.01 


0.84 0.05 


Procyon B*' 


10.82 


-2.28 


13.1 





13.1 0.03 


3.889 0.01 


-1.91 0.01 


0.55 0.02 


Sirius B'^ 


8.44 


-2.89 


11.33 


-2.3 


9.03 0.1 


4.394 0.01 


-2.06 0.01*^ 


1.02 0.02 


40 Eri B° 


9.5 


-1.49 


10.99 


-1.5 


9.49 0.1 


4.223 0.01 


-1.87 0.006 


0.501 0.011 



^from this paper. 

'^from [Harris et al.l (120071 ). 



"from Benedict et all (120001 ). 

^from lO'Brien et all J200lh . BC from lPloweil Jl996l ) and lBereeron et al.l Jl995h . 



''from compilation of jProvencal et al.l (120021 ). BC from iFlowerl (119961 ) and iBergeron et al.l (Il995l ). 
^includes LKH bias correction, negligible for the last four objects. 



sfrom this paper a nd jProvencal et al.l (120021 ), except Feig e 24 from iKawka et al.l (120081 ). Procyon 



B and Sirius B from Schaefer et al. 



ept 

(120061 ) andlBonJ (120091 ). 



^from iBarstow et al.l (120051 ). 
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Table 8. Companion Limits from FGS Fringe Scanning 



ID 


Mv 




M*(©) 


AV 


Comp 2 SpT 


Comp 2M{Q) 


Sep [AU] 


P[y] 


DeHt5 


7.39 


2.9 


0.57 


1 


K8V 


0.59 


3.4 


6 










2 


MIV 


0.53 


5.2 


11 










3 


M2V 


0.4 


17.2 


73 


Abell31 


6.23 


1.61 


0.53 


1 


K4V 


0.7 


6.2 


14 










2 


K8V 


0.59 


9.3 


27 










3 


MIV 


0.53 


31.1 


168 


NGC 7293 


6.77 


4.66 


0.6 


1 


K6V 


0.64 


2.1 


3 










2 


MOV 


0.51 


3.2 


5 










3 


MIV 


0.53 


10.7 


33 


NGC6853 


5.62 


2.47 


0.57 


1 


K3V 


0.72 


4.0 


7 










2 


K6V 


0.64 


6.1 


14 










3 


MOV 


0.51 


20.2 


88 
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Fig. 1. — DeHtS central star and astrometric reference stars. Labels are immediately to the 
right of each star. 
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Fig. 2. — Photometric variations of the four PNNi observed with the FGS. In each panel the 
bottom section shows a fit to the original variation ascribed to roll-induced one year period 
modulation of the stars used to generate each flat field. The top panel shows the residuals, 
each labeled with the final (presumed intrinsic) rms variation in the PNNi. 
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Fig. 3. — Fringes along X and Y axes of DeHtS (solid) and Abell31 (dashed) compared. 
Residuals (bottom traces in each panel) indicate only the amplitude of noise expected for 
these relatively faint targets. 



- 38 - 



I I I I I I I I 



In 
+ 

II 

O 



o 

'a 

o 

T3 



-4 - 



-2 - 



- 



2 - 




0.4 0.8 

J-K (2Mass) 



Fig. 4. — Reduced proper motion diagram for 458 stars in a 1° field centered on DeHt 5. 
Star identifications are shown for DeHt 5 (=1) and our astrometric reference stars. Hk for 
these stars is calculated using our final proper motions (Table Hj). For a given spectral type 
giants and sub-giants have more negative Hx values and are redder than dwarfs in J-K. 
Reference stars ref-9 and ref-11 are confirmed to be giant stars. The cross in the lower left 
corner indicates representative errors along each axis. 
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Fig. 5. — J-K vs V-K color-color diagram. The dashed line is the locus of dwarf (luminosity 
class V) stars of various spectral types; the dot-dashed line is for giants (luminosity class 
III). The reddening vector indicates Ay=1.0 for the plotted color system. 
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Fig. 6. — Histograms of x and y residuals obtained from modeling DeHt 5 and the astro- 
metric reference stars with equations 4 and 5. Distributions are fit with Gaussians whose 
it's are noted in the plots. 
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Fig. 7. — Absolute bolometric magnitude, M^oi plotted against log temperature. The slope 
is constrained to the theoretical value. Mboi for the PNNi comes from our parallax, the 
app arent magnitude, th e bolometric correction, and , the interstellar absor ption, Ay (Table [6]) 
and 



Harris et 



al. 



fjgnnTI ). The WD values are from Benedict et al.l (120001 ) for Feige 24, from 



O'Brien et al.l (l200ll ) for V471 Ta u, and calculated frona Provencal (2002), with bolometric 
corrections from Flower (1996) and lBergeron et al.l (119951 ) for Sirius B, 40 Eri B, and Procyon 
B. All relevant quantities are collected in Table [7l With much scatter the PNNi and WD 
appear to exhibit an approximate absolute bolometric magnitude-temperature relation, one 
that would hold, assuming similar radii and masses for all PNNi and WD. Both Equation 6 
and the residuals, AMhoi, argue otherwise. 
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Fig. 8. — Residuals, AMhoi, from the linear M^oi - log T relationship for PNNi and WD in 
Figure [7] plotted against the logarithm of the radius in solar units. The slope is constrained to 
the value expected from equation 6. The symbol size is proportional to surface temperature. 

The expected linear residual correlation ( ) with log R exhibits an RMS scatter of 0.12 

magnitudes. 
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Fig. 9. — PNNi and WD absolute magnitude, My plotted against log temperature. My for 
the PNNi comes from our parallax, the apparent magnitude, and the inter s tellar absorption, 
Av (Table ED wi th two additional fHD W4, Sh2-216) from Iflarris et al.l J2007h The WD 
values are from [Benedict et al.l (120001 ) for Feige 24, from lO'Brien et al.l (120011 ) for V471 
Tau, and calculated from Provencal (2002). All relevant qu antities are collected in T able [71 
The higher-mass (solid line) evolutionary tracks are fr o m ISchonberner fc Blockerl (119961 ). 
Lower- mass (dashed line) tracks are from iDriebe et al.l (119991 ). The PNNi clump around 
M=0.57Mq. 
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Fig. 10. — PNNi and WD radii from Table [7] plotted agains t log temperature. The evolu- 
tion ary tracks for objects of high and low masses are from ISchonberner fc Blocker! ( 1l996l ) 
and iDriebe et al.l (1l999l ). In this mapping the WD Procyon B and 40 Eri B appear to have 
masses consistent with previous measurements. 
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Fig. 11. — Mass- Radius diagram for six PNNi and five WD with precise radii and masses 
(Table ED. From top to bottom we plot NGC6853, Abell 31, Sh2-216, NGC 7293, DeHt 5, 
Feige 24, 40 Eri B, HDW 4, Procyon B, V471 Tau, and Sirius B. Symbol size is proportional 
to surface temperature. The dashed line is the Hamada & Salpeter (1961) carbon core 
relationship. Except for the coolest PNN, HDW 4, the PNNi (and the hottest WD, Feige 
24) have larger radii than the cooler WD. 
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Fig. 12. — Distances from weighted averages of HST (this paper) and USNO (Harris et 
al. 2007) parallaxes compared to distances from NLTE analysis (wNapiwotzki 2001), H/3- 
derived distances from Shklovski (o 1956), and distances estimated from interstellar Na D 
(A Napiwotzki & Schonberner 1995). The dashed line represents perfect agreement. Objects 
are labeled to the right or top. 



